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ABSTRACT 

A method for computing transient thermal stress vectors from 
temperature vectors is described. The three step procedure 
involves the use of NASTRAN to generate an influence coefficient 
maSx which relates temperatures to stresses in the structural 
model The transient thermal stresses are then recovered and 
sorted for maximum and minimum values. Verification data fo 
the procedure is also provided. 


- 0 INTRODUCTION 

There are many occasions when the stresses produced by transient 
thermal events must be considered. The ascent, on-orbit, and 
descent phases of a Spacelab mission produce large temperature 
aradientson the Cargo Element. A method for recovering the 
thermal stresses produced by these events was developed by the 
Structural Analysis Group at McDonnell Douglas Space Systems 
Company-Huntsvi 1 le (MDSSC-HSV) , and has been used for more than 
six year's in Spacelab Evaluation. 

Because this method was somewhat unhandy to use, only * e^ere*^ 
number of temperature distributions could be run. These were 
generally chosen on the basis of temperature or temperature 
difference extremes. Unfortunately, there is no proven 
relationshipbetween "worst stresses" and "worst temperatures" 
or "worst temperature differences" for the complex m ° d f^ S 

which Spacelab has to deal. It was V ecto?s 

* nrnress that would transform transient temperature vectors 
1° separately calculated) directly into stress vectors 
that could then be sorted for maximum and minimum (Max/Min) 
values . 

Section 2 presents the theory used in the procedure Section 3 
describes the procedure, Section 4 presents verification 
results, and Section 5 contains the conclusions. 
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2.0 


THEORY 


The method described assumes that the thermal strains, and hence 
stresses are linearly proportional to the temperatures. The 
method, therefore, is not applicable to models having 
temperature or stress dependent material properties. 

An influence coefficient matrix, [TRAS], to transform the 
temperatures at a models n grids into m stresses can be 
developed as follows. A temperature of one unit above the 
rsf Grsnce is applied to a grid in the model while the remaining 
n-1 grids are held at the reference temperature. The resulting 
m stresses are then recovered. This is repeated for all n grids 
in the model. The resulting n sets of stresses are then 
assembled as columns to form an m by n influence coefficient 
matrix, [TRAS], that can be used to transform temperature 
vectors into stress vectors. The transformation is as follows: 

[ STMHST] = [TRAS] [TEMPS] (1) 

Where ; 

[STMHST] = Stresses in the finite element model 
(time history) . 

[TRAS] = Linear transformation (influence 
coefficient) from temperatures to 
stresses. 

[TEMPS] = Temperatures at the grids (time 
history) . 

Equation (1) is used to recover transient stress vectors from 
transient temperature vectors. 


3 . o PROCEDURE 

The procedure used to recover maximum and minimum thermal 
stresses is divided into three steps. The first step is to 
generate the influence coefficient matrix [TRAS]. The transient 
stresses are recovered using equation (1) in the second step. 

The transient thermal stresses are then sorted for Max/Min data 
in the third step. Each of these three steps will be described 
below. 
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3 . 1 


Step 1 - Generation of The Influence 
Coefficient Matrix [TRAS] 

NASTRAN is used to generate [TRAS] as described in Section 2. 
NASTRAN Subcase commands are used to accomplish this step with 
each subcase corresponding to a grid in the model. It should be 
noted that NASTRAN can handle a maximum of sixty-six temperature 
load cases. A model having more than sixty-six grids will, 
therefore, require multiple runs and subsequent merging of the 
output data. An example of a NASTRAN control deck to generate 
[TRAS] for the verification model described in Section 4 is 
shown in Figure 1. The DUMM0D5 module is used to convert the 
0ES1 table into a matrix data block which is then written to a 
file using the OUTPUTS module. The extraneous rows (fiber 
distances, safety margins, etc.) of [TRAS] are then removed 
using a specially developed FORTRAN code. 


3.2 Step 2 - Recover Transient Thermal Stress 

Vectors From Temperature Vectors 

The transient temperature vectors must now be obtained. At 
MDSSC-Huntsville, the thermal analyses are performed using 
SINDA. A FORTRAN code has been written to access the SINDA 
output file and recover the desired temperatures along with the 
corresponding time vectors. The temperature and time data is 
written in 0UTPUT5 format. Care should be taken to ensure that 
the row order of the temperature vectors is compatible with the 
column order of [TRAS] . 

Equation (1) is then used to recover the transient thermal 
stress vectors. This is easily accomplished in a simple DMAP 
run using the MPYAD module. The resulting thermal stress 
vectors and time vector are written to a file using the OUTPUTS 
module . 


3.3 Step 3 - Sort Thermal Stress Vectors For 

Maximum and Minimum Data. 

The transient thermal stress vectors are now sorted for Max/Min 
data. This step is performed using a specialized FORTRAN code. 
This code can search multiple time histories allowing composite 
Max/Min tables to be obtained. An example of output from this 
program for the Spacelab Multipurpose Experiment Support 
Structure (MPESS) is shown in Table 1. The data in Table 1 was 
obtained from actual temperature vectors for a Spacelab mission 
and encompasses ascent, orbit, and descent. 
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4.0 


VERIFICATION 


In order to verify that the procedure is working correctly, a 
test case was performed. A simple rectangular plate model was 
constructed using QUAD4 and BAR elements. A plot of the plate 
model is presented in Figure 2 and the Bulk Data is shown in 
Figure 3. The model has thirty-three grids and is homogeneous. 

A temperature gradient (see Figure 4) was applied to the model 
and the resulting stresses in four elements were recovered using 
NASTRAN directly. The results are presented in Figure 5. The 
transient thermal stress recovery procedure was then used to 
calculate stresses due to the same temperature gradient and the 
results are shown in Table 2 . 

It can be seen from Figure 5 and Table 2 that the results are 
the same. This indicates that using [TRAS] to perform the 
thermal stress recovery produces the same results as using 
NASTRAN directly. 


5.0 CONCLUSION 

A procedure has been developed to recover thermal stresses in a 
NASTRAN model directly from temperatures output by a thermal 
model. Because the procedure uses a linear transformation 
matrix rather than a computer program, entire thermal stress 
time histories may be efficiently obtained and scanned for 
Max/Min thermal stress data. Tabular output of the Max/Min data 
is then produced. A test case has been executed and the results 
indicate that the procedure functions correctly. 
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NASTRAN TITLEOPT=0 

s RUN TO GENERATE INFLUENCE COEFF. MATRIX FOR THERMAL STRESS RECOVERY. 

5 

$ 

$ 

$ 

ID CHECKOUT PLATE 
APP DISP 
SOL 1,7 
TIME 999 

DIAC 13,14,21,22,26,42 

5 

$ WRITE OUT ELEMENT STRESS/FORCE MATRICES 
S WHERE; 

$ OESl=STRESSES 

S I 1 OF^ELEMENTS FOR WHICH STRESSES/FORCES ARE BEING RECOVERED 

$ 

DUMMOD5° 6 OESl, , , ,/ELSEU, , , , /C , N , 44/////C , N , l/C.N.l $ 

OUTPUT5 ELSEU, , , , //0/15//0 $ 

OUTPUT5 ,,,,//“ 9/15//0 5 
ENDALTER 
CEND 
$ 

TTTIF = ANALYSIS OF PLATE 

SUBTITLE = RUN TO GENERATE INFLUENCE COEFFICIENT MATRIX 

LABEL = PLATE 
MAXLINES = 99999999 
SPC=100 

$ 

SU LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 1 
TEMP (LOAD) = 1 
ELSTRESS=ALL 

$ 

SU LABEL - APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 2 
TEMP (LOAD) = 2 
ELSTRESS=ALL 

§ 

SU ^ABU - APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 3 
TEMP (LOAD) = 3 
ELSTRESS=ALL 

$ 

SU LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 4 
TEMP (LOAD) = 4 
ELSTRESS=ALL 

$ 

SU LABEL “ APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 5 
Figure 1. NASTRAN Control Deck to Generate [TRAS] or 
Verification Plate Model 
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TEMP (LOAD) - 5 
ELSTRESS-ALL 


- APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 6 
TEMP (LOAD) - 6 
ELSTRESS-ALL 

$ 

SU LABEL - APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 7 
TEMP (LOAD) - 7 
ELSTRESS-ALL 

$ 

SU LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 8 
TEMP (LOAD) = 8 
ELSTRESS-ALL 

$ 

“ ! ARPEV A 1 DEGREE ABOVE REF. TEMP. TO GRID 9 
TEMP (LOAD) = 9 
ELSTRESS-ALL 


SU LABEL =°APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 10 
TEMP (LOAD) = 10 
ELSTRESS-ALL 

$ 

SU L^EL “APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 11 
TEMP (LOAD) = 11 
ELSTRESS-ALL 

$ 

^ LABEL ^APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 12 
TEMP (LOAD) - 12 
ELSTRESS-ALL 

$ 

SU LABEL ^APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 13 
TEMP (LOAD) - 13 
ELSTRESS-ALL 

$ 

^ LAB EL -^APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 14 
TEMP (LOAD) - 14 
ELSTRESS-ALL 

$ 

SU LABEL = 5 APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 15 
TEMP (LOAD) - 15 
ELSTRESS-ALL 

$ 

SU LABEL - 6 APPLY A 1 DEGREE ABOVE REF. for the 

figure 1. NASTRAH^ontrol^e ^ (Continue d) 
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TEMP (LOAD) = 16 
ELSTRESS=ALL 


SUBCASE 17 „ n 

LABEL - APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 17 

TEMP (LOAD) = 17 
ELSTRESS=ALL 

$ 

SU LABEL ^ APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 18 
TEMP (LOAD) = 18 
ELSTRESS=ALL 

$ 

SUBCASE 19 „ rnTn 

LABEL - APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 19 

TEMP (LOAD) = 19 
ELSTRESS=ALL 

$ 

SU !ab!l =°APPLY a 1 DEGREE ABOVE REF. TEMP. TO GRID 20 
TEMP (LOAD) = 20 
ELSTRESS=ALL 

$ 

SU LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 21 
TEMP (LOAD) = 21 
ELSTRESS=ALL 

$ 

SUBCASE 22 ^ „ nT _ 

LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 22 

TEMP (LOAD) = 22 
ELSTRESS=ALL 

$ 

SUBCASE 23 „ nTT _ 0 - 

LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 23 

TEMP (LOAD) =23 
ELSTRESS-ALL 

$ 

SU LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 24 
TEMP (LOAD) = 24 
ELSTRESS=ALL 

$ 

SUBCASE 25 „„ „ T _ 

LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 25 

TEMP (LOAD) “ 25 
ELSTRESS=ALL 

$ 

SUBCASE 26 

LABEL = APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 26 

TEMP (LOAD) = 26 

ELSTRESS-ALL 

$ 

SUBCASE 27 

LABEL - APPLY A 1 DEGREE ABOVE REF. TEMP. TO GRID 27 

Figure X. NASTRAN Control Deck to Generate [TRAS] fo 
Verification Plate Model (Continued) 
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TEMP (LOAD) « 27 
ELSTRESS'ALL 


SUBCASE 28 

LABEL = APPLY A 
TEMP (LOAD) = 28 
ELSTRESS=ALL 

$ 

SUBCASE 29 

LABEL = APPLY A 
TEMP (LOAD) = 29 
ELSTRESS=ALL 

$ 

SUBCASE 30 

LABEL = APPLY A 
TEMP (LOAD) = 30 
ELSTRESS=ALL 

$ 

SUBCASE 31 

LABEL = APPLY A 
TEMP (LOAD) = 31 
ELSTRESS=ALL 

$ 

SUBCASE 32 

LABEL = APPLY A 
TEMP (LOAD) =32 
ELSTRESS=ALL 

$ 

SUBCASE 33 

LABEL = APPLY A 
TEMP (LOAD) =33 
ELSTRESS=ALL 

$ 


1 DEGREE ABOVE REF. TEMP. TO GRID 28 


1 DEGREE ABOVE REF. TEMP. TO GRID 29 


1 DEGREE ABOVE REF. TEMP. TO GRID 30 


1 DEGREE ABOVE REF. TEMP. TO GRID 31 


1 DEGREE ABOVE REF. TEMP. TO GRID 32 


1 DEGREE ABOVE REF. TEMP. TO GRID 33 
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Figure 1. NASTRAN Control Deck to Generate [TRAS] for the 
Verification Plate Model (Concluded) 
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Figure 3. Plate Model Bulk Data (Concluded) 
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Figure 4 . Temperatures Applied to the Plate Model 
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Figure 5. Thermal Stresses from a Direct NASTRAN Solution of 
the Plate Model 
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Figure 6. MPESS Finite Element Model 
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Figure 7. MPESS BAR 106 (SINDA Node 6000) Axial Stress Time 
History 
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Figure 8. MPESS BAR 106 (SINDA Node 6000) Temperature Time 
History 
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MDSSC - HUNTSVILLE OPERATIONS 


TRANSIENT THERMAL STRESS RECOVERY PROCEDURE 
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